Black Hole in the Expanding Universe with Arbitrary Power-Law Expansion 
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We present a time-dependent and spatially inhomogeneous solution that interpolates the extremal 
Reissner-Nordstrom (RN) black hole and the Friedmann-Lemaitre-Robertson- Walker (FLRW) uni- 
verse with arbitrary power-law expansion. It is an exact solution of the D-dimensional Einstein- 
"Maxwell" -dilaton system, where two Abelian gauge fields couple to the dilaton with different cou- 
pling constants, and the dilaton field has a Liouville-type exponential potential. It is shown that 
the system satisfies the weak energy condition. The solution involves two harmonic functions on 
a (D — l)-dimensional Ricci-flat base space. In the case where the harmonics have a single-point 
source on the Euclidean space, we find that the spacetime describes a spherically symmetric charged 
black hole in the FLRW universe, which is characterized by three parameters: the steepness pa- 
rameter of the dilaton potential Ut, the U(l) charge Q, and the "nonextremality" r. In contrast 
with the extremal RN solution, the spacetime admits a nondegenerate Killing horizon unless these 
parameters are finely tuned. The global spacetime structures are discussed in detail. 

PACS numbers: 04.70.Bw,04.20.Dw,04.50.+h, 04.50.Gh 
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I. INTRODUCTION 

Black holes could have formed in the early stage of the 
universe as a consequence of primordial density fluctu- 
ation P, 0] ■ Since these primordial black holes (PBHs) 
are of the order of the horizon scale mass, they span 
a considerably broad mass spectrum ranging from the 
Planck mass (10 _5 g) to 1O 5 M0, or even much larger. 
Hence, their presence could leave diverse physical im- 
prints throughout the cosmic history (see e.g., Q and 
references therein). The PBHs with M > 10 15 g survive 
until the present epoch, so that they are plausible can- 
didates of cold dark matter, and likely origin of super- 
massive black holes and/or the sources of gravitational 
waves. The PBHs with mass M ~ 10 15 g are now evap- 
orating via Hawking evaporation Such PBHs emit 
quanta of order lOOMeV, which could also contribute to 
the cosmological 7-ray background and generate 7-ray 
bursts. While, the PBHs with mass smaller than 10 15 g 
that have completely evaporated during a first second of 
big bang could generate large amount of entropy of the 
universe. Thus, the number density of small PBHs could 
place strong constraints on the big-bang nucleosynthesis. 
In these contexts, studying the formation and the evap- 
oration process of black holes in the expanding universe 
are of great significance as a probe of the early universe, 
high energy physics, and quantum phase of gravity. 

Apart from the above astrophysical interests, a study 
of black holes in the expanding universe has been a prin- 
cipal stirring subject in general relativity. Black holes 
have played a central role in general relativity, since they 
encode the essential characteristics of the gravity theory, 
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viz, the nonlinearity and highly curved spacetime. Vari- 
ous kinds of physical and geometrical properties of black 
holes have been clarified thus far. The most enduring 
achievements among these is the uniqueness theorem of 
black holes 0, @ , according to which isolated black holes 
in equilibrium states-these states occur if sufficiently long 
time passed after the gravitational collapse of massive 
stars-necessarily belong to the Kerr family. This means 
that stationary black holes are completely characterized 
by two conserved charges (mass and angular momentum) 
without any additional "hairs." By virtue of this theo- 
rem, stellar-sized black holes with futile ambient sources 
are well approximated by Kerr black holes. To the con- 
trary, in the non-isolated and dynamical background, a 
variety of black holes equipped with much richer proper- 
ties are expected. However, when one attempts to obtain 
the black hole spacetime in the dynamical background, 
a serious difficulty arises. If a black hole is put on the 
homogeneous and isotropic FLRW universe on which the 
standard cosmological scenario lays the foundation, the 
background universe will become inhomogeneous, and at 
the same time the black hole will continue to grow and/or 
deform by swallowing ambient matters. Hence, the fact 
that the spacetime is time-evolving and spatially inho- 
mogeneous enforces us to solve nonlinear partial differen- 
tial equations for the geometry, as well as for the matter 
fields. 

The exact black-hole solutions in the FLRW universe 
found hitherto have enjoyed high degrees of symmetry. 
Among them are the Einstein-Straus model [3 (a black 
hole in the "Swiss-Cheese Universe"), the Schwarzschild- 
de Sitter (SdS) and the Reissner-Nordstrom-de Sitter 
(RNdS) black hole. Since these solutions maintain the 
equilibrium due to the timelike Killing field, they are 
unlikely to capture the envisaged dynamical fact in re- 
alistic situations. In recent years, Sultana and Dyer 
have exploited the conformal technique and have con- 



2 



structed a dynamical black hole [8[ which is asymp- 
totically Einstcin-dc Sitter universe and conformally re- 
lated to the Schwarzschild solution. Unfortunately, the 
Sultana-Dyer solution suffers from violating energy con- 
ditions. Meanwhile, if one imposes a self-similarity (char- 
acterized by a homothetic Killing field), it has been 
proved that a spherically symmetric black hole cannot ex- 
ist in the asymptotically decelerating FLRW universe @ . 
These examples illustrate the difficulty in constructing a 
regular black hole immersed in the FLRW universe with 
desirable physical properties. 

Recently, a time-dependent "black hole candidate" was 
obtained via the dimensional reduction of the dynam- 
ically intersecting branes in 11-dimensional (11D) su- 
pergravity The solution appears to behave like a 

charged black hole for small radii, while for large radii it 
approaches to the FLRW universe filled with stiff mat- 
ter fluid. In the previous paper we elucidated the 
detailed spacetime structure and established that the so- 
lution indeed describes a charged black hole that ap- 
proaches to the flat FLRW universe. It is shown that 
the metric is an exact solution in Einstein- "Maxwell" - 
dilaton system with four kinds of U(l)-ficlds (three of 
them are degenerate) coupled to the dilaton, thereby 
the system satisfies the dominant energy condition. Al- 
though the BPS black hole plays the leading part in string 
theory, what is intriguing us is that the solution is not 
extremal, i.e., the Hawking temperature does not van- 
ish. Despite the little astrophysical concern because of 
the brane charges, the intersecting brane picture will en- 
able us to obtain a variety of dynamical black holes with 
physical matter sources. In [12|], Gibbons and one of the 
present authors found the "black hole candidate" that 
asymptotically looks like an FLRW universe with arbi- 
trary power-law expansion by introducing an exponen- 
tial potential. The solution found in [l2[ reduces to the 
solution in [ljl [H| as a special case and is expected to de- 
scribe a black hole in the expanding universe. As lessons 
from the McVittie solution (l3l - [l5| . a great deal of care 
is needed in order to obtain the global spacetime struc- 
ture of a time-evolving and spatially inhomogeneous so- 
lution. 1 In this paper, following the previous study (Tl| 
we intend to clarify the spacetime structure of their so- 
lution, which has been an open issue in [l2j . 

There are several motivations for studying black holes 
in the Einstein- "Maxwell" -dilaton system with an expo- 
nential potential. In the light of string theory, the dila- 
ton and the form fields are elementary constituents in the 
theory, and the exponential potential of a dilaton natu- 
rally arises in various contexts: the excess central charge 
(the "noncritical" string theory) [l7|, the massive IIA 
supergravity [Tc| , the Kaluza-Klein compactification of a 



1 In contrast to the claims in 1 1 -ll , it is recently argued that McVit- 
tie's solution includes a regular black-hole horizon if the positive 
cosmological constant dominates the universe at late time [Tr3l . 



curved internal space [19( , the vacuum expectation value 
of the four-form field strength [2(j , and so on. From the 
general relativistic point of view, these fields obey suit- 
able energy conditions. Thus, it is interesting to see if 
the black hole exhibits a peculiar aspect under a realis- 
tic circumstance. In addition, since the exponential po- 
tential may drive the power- law inflation [21[ , this could 
give a profound implication for the background geometry 
of PBHs during inflation. Furthermore, by generalizing 
the present solution to include multiple black holes one 
is able to discuss the collision and coalescence of black 
holes, providing a valuable arena to test the cosmic cen- 
sorship conjecture (22|. 

This paper constitutes as follows. In the ensuing sec- 
tion, we present a D-dimensional solution as a simple 
and honest generalization of the solution given in [1 21 ] . 
We discuss the matter fields and singularities of the so- 
lution in Section IIII1 where the system is shown to sat- 
isfy the suitable energy conditions. Section ITVl develops 
our main argument on the spacetime structure. Analysis 
of the near-horizon geometry and null geodesic motions 
leads to the conclusion that the present spacetime can 
admit a regular event horizon with constant circumfer- 
ence radius, if the parameters of the solution are chosen 
appropriately. Conclusions and future outlooks are sum- 
marized in Section fVl 

We shall work in units c = K = 1 with k 2 = 8ttG and 
follow the standard curvature conventions 1Z? vpa V v := 
2V [p V (T] V^, := Wppv and K := Since the 

circumference radius will be denoted by R, the script no- 
tation is used for curvature tensors throughout the paper. 



II. PRELIMINARIES 

The solution describing a charged black hole in the 
FLRW universe is derived in fioj from the dimensional 
reduction of dynamically intersecting M2/M2/M5/M5- 
branes in I ID supergravity. Viewed from 4D inhabi- 
tants, the solution satisfies the field equations in Einstein- 
"Maxwell" -dilaton system with particular couplings. The 
gauge fields and the scalar field arise via toroidal com- 
pactification of extended directions of branes in 11D su- 
pergravity. The massless scalar field is responsible for the 
stiff-fluid dominant universe. This solution was general- 
ized in [l2j in such a way that the background FLRW 
cosmology obeys the arbitrary power-law expansion by 
introducing the Liouville-type potential. 

In this paper, a more general class of solutions are 
to be considered. We shall begin by the TD-dimensional 
Einstein- "Maxwell" -dilaton system, in which two types 
of U(l)-fields couple to the dilaton with different cou- 
plings, and the dilaton has a Liouville-type exponential 
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potential. To be specific, the action is described by 
1 



S 



d D x- 



-fJ 



2k 



1D7T ^-^ 



1 

2 



with 



= V r exp(-«K$) . 



(2.1) 



(2.2) 



Here, a (> 0) is a dimensionless constant corresponding 
to the steepness of the potential. We have introduced 
degeneracy factors, ua (> 0), of two U(l) fields for later 
convenience (which may be absorbed into the definition 

of F^). This type of action with a single U(l) field 
has been discussed in some situation, e.g. a black hole 
collision in the D = 4 theory with a = — At = 2 [23[ (see 
also [13). 

The above action (|2.1[) yields the field equations for 
the metric, the dilaton, and two U(l) fields as, 



K 



dV 
F d<I> 



„-^ V = « 2 (T«+jJ»)) , (2.3) 

. l - V n A X A Ke x -^F^F^ = , 
Ion p 

(2.4) 

where 

- (V„*) (V„$) - ^ [(V*) (V p $) + 2V] , 

(2.6) 
(2.7) 



nn(em) 



E 



4tt 



Inspired by the solutions given in [10|, [l2| , let us specify 
the parameters (tia, \a,o:) as 



n-T + n s 
Xt = a = 

As = -2 



2{D-2) 
D-3 



'(J-3)n s 
(D - 2)n T 



(£)-3)n T 



(2.8) 
(2.9) 

(2.10) 



Y (£> - 2)n s ' 
for which 

0<77, T <4, 0<n<?<4, A T >0, \ s < . (2.11) 

The constants tit and ns may take natural number only 
for D = 4, 5, in which case they are related to the num- 
ber of time-dependent and static branes in 11D super- 
gravity (Tol[Tl|. In the present case, however, tit and ns 



are not necessarily integers, hence they do not have such 
a meaning. 

With the above choice of parameters, we shall look for 
a spatially-inhomogeneous and time-dependent solution 
of the following form, 



da* 



-E D - 3 dt 2 + S- 1 da 2 D _ x , 



with 



(1 + H S ) 

hij (x)dx I dx J , 



-1/(73-2) 



(2.12) 

(2.13) 
(2.14) 



where to is a constant with dimension of time, and Ht(x) 
and Hs{x) are functions of x 1 on the (D — l)-dimensional 
base space da 2 



D-V 



Assuming that the base space is Ricci-flat, 

h n u = o, 



(2.15) 



where IZij is the Ricci curvature of the spatial metric 
hu, and Ht{x) and Hs(x) are harmonic functions of the 
base space da 2 D _ 1 , 



h AH T = , h AH s = , 
and if t is tuned to be 



nxin-T — 1) 
4k 2 Vq 



(2.16) 



(2.17) 



we find that the metric (|2.12|) solves the field equations 
(|2.3[) - (|2.5[) provided that the dilaton and the electromag- 



netic potentials defined by 



are given by 



1 (D-3)n T n s 



D-2 



In 



kA[ T) = V^Ht 1 + a (T) (t), 



kA 



(S) 



2nH s L +a^(t) , 



where 



Ht — 



— + Ht , Hs — 1 
to 



Hs 



(2.18) 



(2.19) 

(2.20) 
(2.21) 



(2.22) 



and a^ T '(t) and a^ s '(t) are arbitrary functions of t. 

Bearing the cosmological application in mind, we set 
hu to be a fiat Euclidean metric, hu = Sjj, in what 
follows. Note that the action is not invariant under elec- 
tromagnetic duality transformation except the = 1 
case. We shall thence focus on the electrically charged 
solution throughout the article. 
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Let us consider the case where the harmonics represent 
a system of N point-sources, 



N n (i) 



i=l 
N 



x —x 



U) 



(2.23) 
(2.24) 



where and Qg are charges at the points xu\ (i = 
1, • • • ,N), respectively. It follows that the case of 
?it = reduces to the higher-dimensional Majumbdar- 
Papapetrou solution J2^|, while the case of 715 = 
describes the higher-dimensional Kastor-Traschen solu- 
tion (2(| [27j • Henceforth, we do not consider these spe- 
cial cases unless otherwise stated, i.e., we assume nx 7^ 
and ns 7^ 0. For D = 4, tit = 1, one recovers the solution 
in [l(| , which has been shown to describe a charged black 
hole in the FLRW universe when each harmonic function 
has a single point source at the origin (Tlj . 

Supposed (t/to) > 0, let us transform to the new time 
coordinate t defined by 



t 



(P-3) 

2TD^27™S _ 2(D-2) 

with t = 7^ ./ t , 2.25) 



in terms of which we can cast the metric (|2.12[) into the 
form 



ds 2 = -Z D - 3 dP + a 2 E- 1 S IJ dx I dx J 



(2.26) 



where 



1 2(X)-2)/n T 



(i+H S y 



-l/(£»-2) 



(2.27) 



and 



with p = 



Tlx 



(D - 3)n s 



(2.28) 



Since we imposed the boundary condition such that the 
harmonics Ht and Hs fall off as r := \/J2i( xI ) 2 ~^ 
00 [Eq. (|2~23|) and (|2~24)) ]. the metric ([2~26]) approaches 
in the limit r — > 00 to the £>-dimensional flat FLRW 
spacetime, 



d.s. 



-dt 2 + a 2 5 7J da; 7 da; J 



(2.29) 



The new coordinate i is found to measure the proper time 
at infinity. Looking at the behavior of the scale factor 
a oc f p , one can recognize that the asymptotic region of 
the spacetime is the FLRW universe filled by a fluid with 
the equation of state 



P 



top, 



with 



2(D - 3)n s 



- 1 



(2.30) 



(D-l)n T 

It turns out that the parameter n<r (or ns) is associated 
to the expansion law of the universe (|2.28[) . Notably, we 
can obtain an accelerating universe (p > 1) by setting 
riT > 2 or equivalently ns < 2/(£) — 3). In particular, 
the exponential expansion (the de Sitter universe) is un- 
derstood to be p — > 00 (ns —> 0). Figure [1] depicts the 
conformal diagrams of the FLRW universe. The asymp- 
totic regions of the present spacetime (|2.12l) resemble the 
corresponding shaded regions in Figure [T] 




FIG. 1: Conformal diagrams of a fiat FLRW universe a = (t/t ) p for (1) < p < 1/2, (2) p = 1/2, (3) 1/2 < p < 1, 
(4) p = 1 and (5) p > 1. The dotted and dotted-dashed lines denote the trapping horizon, rTn{t) = (da/df) -1 , and the 
big-bang singularity at a — 0, respectively. The cases (2) and (4) correspond respectively to the radiation-dominant universe 
P — p/(D — 1) and the marginally accelerating universe driven by the curvature term p oc a~ 2 . The cosmological horizon, 
?*CH(t) = {p — l) -1 io(i/fo) 1 ~ p , is abbreviated to CH, which exists only in the strictly accelerating case (p > 1). The shaded 
regions corresponding to r — > 00 approximate our original spacetime. 
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r = 0,t= +<*> 
t= const. <oo 




r=0, t= const. 



FIG. 2: Typical "throat" geometries for the AdS 2 X S D ~ 2 
spacetime (left) and the extremal RN spacetime (right). The 
double line stands for the degenerate horizon. Only the vicin- 
ity of the shaded point at r — with t kept finite approximates 
our original solution. Whereas, the domains bounded by blue 
shaded lines in each figure are locally isometric [f|. 



On the other hand, taking the limit rt^ := \x — xu)\ — > 
0, we can safely neglect the time-dependence of the met- 
ric. Hence, the metric at the very neighborhood of 
each mass point is approximated by the Nariai-Bcrtotti- 
Robinson metric [HI, [29[ [the direct product of a 2- 
dimcnsional anti-de Sitter (AdS2) and a (D — 2)-spherc] 
as 



ds 



pi al 



ft 

'All 

r 2 



dr 



(0 



D-2 > 



(2.31) 



where £ {l) := [(Q^ ) ) n ^(Q^ ) ) n s]V[2(c-2)] setg the curva _ 

ture scale of AdS2 and S 2 at the i-th point, and dfl 2 J _ 2 is 
the line-clement of a (D — 2)-dimcnsional unit sphere. It 
has been noticed that the neighborhood of any extremal 
black holes can be universally described by the above 
metric [10, HH . Figure [5] compares the geometry of the 
AdS2 x S D ~ 2 and that of the extremal RN black hole. 

Thus, one anticipates that the metric (|2. 12|) describes 
a system of charged black holes with a degenerate event 
horizon embedded in the FLRW universe filled by a 
fluid (|2.30p . which might lead us to speculate that the 
causal structure is found simply by patching Figures Q] 
and [2] As pointed out in [ll|, however, this naive expec- 
tation fails if a single point mass is placed at r = for 
D = 4 and ut = 1. The solution in [k| has neither a 
degenerate event horizon at r = 0, nor a big-bang sin- 
gularity at t = 0. As is clear from Figure [2J the limit 
r — > with t kept finite does not describe the (D — 1)- 
dimensional null surface (the blue shaded domain corre- 
sponding to r — > with \t\ going to infinity), but de- 
scribes only the point at the "throat infinity." This lies 
at the core of why the present metric admits nonvanish- 
ing Hawking temperature. In order to elucidate the full 



causal structure, one has to discuss these null surfaces 
with much care, the analysis of which is the principal 
subject of this paper. 

A curious property of the present solution (|2.12[) is 
that the field equations are completely linearized. One 
may encounter the linear gravitational equations in the 
context of the solitonic configurations and the supersym- 
metric solutions. In the present case, none of these inter- 
pretation fails. Despite the non-trivial time dependence, 
we are able to superpose the harmonics ad arbitrium. 
This property might be seemingly contradictory. Yet, 
one faces the same situation in the Kastor-Traschen solu- 
tion [2^] , where the black hole collision can be described 
by multicenter harmonics. These spacetimes maintain 
a special balance in close analogy with supersymmetric 
states. 2 

We will postpone the discussion regarding the multiple 
black holes in a separated paper. In this article, we shall 
limit our attention to the simplest case in which the har- 
monics Has possess only a monopole term cx QaIt DZ 
at the coordinate center, i.e., the spacetime is spherically 
symmetric. In this single-mass case for each harmonic, 
it is sufficient to consider the case of identical charges 
Qt = Qs ='■ Q; since Qt 7^ Qs amounts to a trivial 
conforms! change 



9^ (Qt/QsT t,(D - 2) 9^, 
with the parameter redefinitions 

* -> (Q T /Qs)- nT/2 t, 
to^ {QT/Qs) l - nT,2 t . 

Thus, we will be concerned with the metric 
ds 2 = -E D - 3 dt 2 + S- 1 (dr 2 + r 2 di " 

where 



D-2 ' 



with 



Ht \ — — 



H, : = 1 



r-D-3 



Q 



r.D-3 



(2.32) 

(2.33) 
(2.34) 

(2.35) 

(2.36) 

(2.37) 
(2.38) 



2 It has been pointed out that in the Kastor-Traschen spacetime 
there exists a super-covariantly constant spinor by Wick-rotating 
the cou pling c onstant, g, in the gauged N = 2 supergravity as 
g — > ig [27ll32j . This is not related to the genuine supersymmetry 
because gauged supergravities require the negative cosmolo gica l 
constant. Rather, it arises from the "fake" supergravity [33-35], 
in which the supercovariant derivative operator is no longer Hcr- 
mitian, and the unwanted ghosts appear due to the non-compact 
gaugings. It is interesting to discuss if the present system can be 
embedded in the fake supergravity. 
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The constant Q has a physical meaning as a U(l) charge, 
satisfying 



where S is an arbitrary closed surface surrounding r = 0, 
and 



&D-2 



2% (D-i)/2 
r[(£»-l)/2] 



(2.40) 



denotes the area of a (D — 2)-dimensional unit sphere. 
The above first integral is a direct consequence of 
Eq. (|2.5p . Though the charge may take either sign, we 
shall restrict ourselves to the positive case (Q > 0) for 
definiteness of our argument. 



III. MATTER FIELDS AND SINGULARITIES 

Let us first evaluate matter fields and discuss energy 
conditions. Next, we shall discuss the spacetime singu- 
larity. 



A. Energy conditions 

Many candidates of "black hole solutions" in the ex- 
panding universe found in the literature (Sultana-Dyer 
solution [H, McVittie's solution [l3j . etc) do not respect 
energy conditions in the whole of spacetime. It is there- 
fore worthwhile to argue whether the present spacetime 
satisfies suitable energy conditions. 

Since the explicit metric-form is available, we can im- 
mediately evaluate the matter contents of the solution. 
Let us consider a (nongcodesic) observer whose integral 
curve is tangent to d/dt. The energy densities and the 
pressures for this observer are given by 



p(*) 



i r 



s-( D - 3 )($ t ) 2 + s($ r ) 



p. 



Ti + T 2 + T 3 > , 
I ^F,-( D - 3 )(<t>A 2 + P, 

l r. 

2 

J~l — J 7 2 — J 7 3 



S -P-3) ($it) 2 + H($ir) : 
T2-T3, 

5-(^-3)($ t )2_ S ($ r) 2 



V 



V 



V 



(em) _ _ p(cm) _ p(em) 



-P, 

T 4 > 0, 



(3.1) 

(3.2) 

(3.3) 
(3.4) 



where the comma denotes the partial derivative, and 

* S (P-3) 



Fa 

^4 



n T n s (D-3) ( H, 



'8/t 2 (D-2)t 2 \H T 



(D-2) 



n T n s Q 2 {D-Zf{H T -H s Y 



t 2( D _ 2 ) r 2(D-2)p, 



^ (D-2) XT (-D-2) 

T n S 



■» S ( p -3) 
(D-2) 



nxiriT — 1) 
4^ 

(D - 3) 2 Q 2 (n T H 2 + n s H 2 ) 
AK 2 r 2(D-2) H 2 + (S) H 2+ A 



(3.5) 
(3.6) 

(3-7) 
(3.8) 



Here, the subscript D, represents any of the angular com- 
ponents, which arc all equal under the spherical symme- 
try. Due to the scalar field, the observer sees the non- 
vanishing energy flux, whose radial component T§ := 
_ E -(D-4)/2 Ttr is found tQ be 



= - 



n T n s Q(D-3) 2 {H™ T H r s 



(D-4) 

l (d-2) (h t _ h s ) 



4k 2 (D - 2)t r D - 2 H 2 H s 



(3.9) 



From these, one can derive the following useful relations 



•Fi + F 3 = 
F 4 - T 2 = 



n s (D-3) 
(D-2) 



n 2 T (D-l) fH 
8n 2 (D-2)t 2 \H T/ 

(D - 3) 2 Q 2 ffg T 

8k 2 (D - 2)r 2 ( D - 2 )H 2 T +T ^ 1T H 2 s +T ^ 1J 



(3.10) 



where we have used a shorthand notation 
Hst : = u sHt + u tHs ■ 



(3.11) 
(3.12) 

(3.13) 



Note that T\, T2, J 7 4, T\ + J-3 and T4 — Ti are positive- 
definite, while J- 3 is negative for < ut < 1. 

At the time t = to we have Ht = Hs, hence the physi- 
cal quantities evaluated at t = to are reminiscent of those 
for the extremal RN solution. We shall refer to the time 
t = to as a fiducial time. At t = to, the energy density of 
the scalar field is uniform in space as 



2 ($) _ ( D - l ) n T 
KP ~8(£-2)t 2 ' 



(3.14) 



and that of the electromagnetic field is the same as the 
extremal RN solution, 



(D-2)(D-3)Q 2 

2( r D-3 + Q)2(£>-2)/(D-3) 



(3.15) 



For t > t , J-5 becomes negative, indicative of the scalar 
field falling into the black hole. As it turns out in the 
later section, this is not the case since the black-hole 
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event horizon is a Killing horizon. Killing horizons rep- 
resent "equilibrium states" of event horizons, precluding 
the black hole to grow [36]. Physically considering, the 
accretion of the scalar field is delicately in balance on 
the horizon with the repulsive forces caused by the U(l) 
fields. 

We are now ready to discuss energy conditions for the 
total energy-momentum tensor. For the observer moving 
along d/dt, the total energy density and the pressures 
are given by 



,(*) 



n( Gm ) 



Pr = + P r (Cm) , 



p< 



(cm) 



(3.16) 
(3.17) 
(3.18) 



The simplest and the most intuitive one is the weak en- 
ergy condition [37| . which states that the energy den- 
sity for an arbitrary observer with worldline is non- 
negative, 



In terms of these quantities, the weak energy condition 
reads (H 



P>0, p + P r >0, P + Pn>0. (3.23) 



From Eqs. (|3.11[) and f|3.18[) one obtains 



[{p + P r f-^Fl] x ' 2 = 2\T X -T 2 



(3.24) 



Noticing J- 4 — J- 2 > 0, it follows that the last two inequal- 
ities in Eq. (|3.23[) hold. Then, it remains to examine the 
positivity of p = F 3 + F A + \Fi — or a lack thereof. 
In the region T\ — J 7 2 > 0, we have 



Tu V V^V" > . 



(3.19) 



Solving the characteristic equation dct(T A '^ — T5^ u ) = 
0, we can diagonalize the (total) stress-energy tensor as 
T^p = diag(— p, P r ,Pn,...,Pn), where the eigenvalues are 
given by 



„ 1 
P= 2 



^=2 

Pa = P 



p-P r + J(p + P r y-4Ti 



(p - P r ) + J (p + P r f - 4Ti 



n ■ 



(3.20) 

(3.21) 
(3.22) 



p = {Fx + J" 3 ) + (^4 - T 2 ) > , 



while in the region J-\ — J- 2 < 0, we have 



p = —T\ + T2 + F3 + ?4 ■ 



(3.25) 



(3.26) 



At first sight, the sign of J-3 + J- \ might be obscure in 
the latter case. Yet, using the condition T\ — J- 2 < 0, it 
turns out to be bounded below as 



J3 + J 7 4 > 



riTHg + nsHx + tit (pit 



1){Ht-H s )^ 



(D-2) 



~{n s 



2nx Hx H s 



n 2 T (H T 



H s f 



(3.27) 
(3.28) 



From Eq. (|3.27l) the positivity of p immediately follows 
for jit > 1. In the case of nr < 1 [ns > (D— 1)/(D — 3)], 
observing the inequality 



D - 1 

n?- nr > 1 = 

6 D-3 



D 



>0, 



(3.29) 



we are led to the conclusion p > from Eq. (|3.28|) . Hence, 
the weak energy condition is satisfied in either case. 

One can evaluate other types of energy conditions in 
a similar fashion. In the background FLRW universe, 
the acceleration of the universe is driven when timclikc 
convergence condition ceases to be valid [37| ■ In general 



relativity this condition is equivalent to saying that the 
strong energy condition 



1 



D-2 



T p P g^ VV > , (3.30) 



is violated, indicative of peculiar matter fields with large 
negative pressures. In terms of eigenvalues of the stress- 
energy tensor, the strong energy condition amounts to 

p + P r > , p + Pn > , 

{D-3)p + P r + (D-2)Pn > 0. (3.31) 
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The last condition is rewritten as 

(D - 3)p + P r + {D- 2)P n 
= {D-2){\F l -F 2 \+F 1 -F 2 ) 
+ 2 [- J- 3 + (D - 3) J4] . (3.32) 



Since the last term on the right-hand side of Eq. Q3.33j) 
is always dominant for r — > 00, the strong energy condi- 
tion violates if ut > 2, i.e., when the background FLRW 
universe is accelerating (p > 1). While, for ut < 2, 
i.e., when the background FLRW universe is decelerating 
(p < !)■ wc can conclude that the strong energy condi- 
tion holds, following the same line of argument as Eqs. 
(|3.25|) — (|3.29|) (note that we must also evaluate the case 
of T\ — T2 < 0). Namely, the present spacetime satisfies 
the strong energy condition if and only if the background 
FLRW universe does. 

As shown above, it is dependent on the parameters and 
restricted to the particular spacetime regions whether the 
strong (and the dominant) energy condition holds. The 
failure of these energy conditions is, however, not so se- 
rious. It is widely accepted that the minimal require- 
ment to avoid any pathologies is the null energy condi- 
tion: p + P r > and p + Pa > 0. If the null energy 
condition is not satisfied, undesirable instabilities would 
develop Since the present system satisfies the weak 
energy condition, the null energy condition automatically 
holds H?! • Hence, the present spacetime is free from such 
annoying plagues. This motivates us to study the system 
as a physically acceptable one. 

B. Misner-Sharp energy 

A celebrated quantity that characterizes the matter 
fields in spherical symmetry is the Misner-Sharp quasilo- 



The present spacetime metric (|2.12[) yields 



In the case of J-\ — J- 2 > 0, a straightforward calculation 
reduces the right-hand side of the above equation to 



(3.33) 

I 

cal energy [3^| . The advantage of the use of Misner-Sharp 
energy is that it also characterizes the local spacetime 
structure (40l - [42l |. The quasilocal energy is interpreted 
as a gravitational energy contained within a closed sur- 
face. Especially when the spacetime fails to admit glob- 
ally conserved charges, the quasilocal energy is a helpful 
local measure to identify matter distributions. 

A major superiority in spherical symmetry lies in the 
fact that we can define covariantly the circumference ra- 
dius as 



R := r~- 1/2 (3.34) 

Using the above circumference radius, the D-dimensional 
Misner-Sharp quasilocal energy, m, is given by (4l| 



m = ^—^Q D - 2 R D ~ 3 [1 - ff""(V M .R)(V„iZ)] . (3.35) 



The Misner-Sharp mass is defined uniquely by the metric 
components and its first derivatives, and does not require 
the asymptotic conditions of spacetime. Once the spher- 
ical surface specified by the circumference radius is fixed, 
the Misner-Sharp mass is given without ambiguity. Such 
a quasi-localization is possible on account of spherical 
symmetry, in which no gravitational wave is generated. 



8t 2 oK 2 H 2 T +( °- 2) H 2 s +( °- 2) 



{D r 2(D-?) ~ { 2 ( n TH 2 s + n s H 2 ) - n T n s (H T - H s ) 2 } - (D - l)n T (n T - 2)Hp H n s s+2 



n D -2\r\ D - 3 

8(D - 2)t2 K 2 r 2(£>-l)^2-[(Z>-3)n T ]/[2(D-2)] ff 2-[(Z3-3)n s ]/[2(D-2)] 

x [n 2 T r 2D H^ T H n s s+2 + 4(D - 2){D - 3)Qr D+1 t 2 H T H s H S T - (D - 3) 2 Q 2 r 4 t 2 H 2 ST ] . (3.36) 



We can provide each term in Eq. (|3.36[) with physical interpretation by evaluating at the fiducial time t = to, at which 
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the Misncr-Sharp energy reduces to 



m (to) R) = l (D _ 2 ^ 4 + ^iT s ^Q + -tar"*-* 



£=i^-V*5(*o) + n D _ a ^ dRR?-*p™ (to, R) + {D 2) f J - 2 Q . (3.37) 

- 1 7 Q i/cr>-3) 2k z 



fi 

D-l 



The first term in the above equation corresponds to the 
energy of the scalar field, and the last two terms to the 
U(l) energies outside and inside the black hole. 

For Q = 0, only the first term in Eq. (|3.36[) survives, 
yielding 



as expected for the background FLRW universe. This 
implies also that the first term in Eq. (|3.36j) is the con- 
tribution of a scalar field. 



(D-2)0 D _ 2 a i 



2k 2 



D-l 



(ar) 1 



-y*) 



(3.38) 



C. Curvature singularities 



We can see easily that the dilaton profile (|2.19j) and 
the square of the electromagnetic fields 



f (t )f{ t)„ = MD-z) 2 Q 2 ( h? 



n s (D-i) \ UlD-2) 



r 2(D-2) I ^(4-n T )-D+4(n T -2) 



(3.39) 



r L LV r 2{D-2) \ H (A-n s )D+4(n s -2) J ' \0-*V) 

diverge at 

t = t s (r) := -Q/r D - 3 , and r°- 3 = -Q . (3.41) 

One may thence deduce that these surfaces correspond to the genuine spacetime singularities characterized by the 
divergence of scalar quantities of curvature. As a simple illustration, let us see the square of the Weyl-tensor 

C^ pff C^ = |^-^W 2 , (3.42) 

where (see e.g., appendix of [4l| for derivation) 

(D - 3)Q [QD - 3)Qr 3 {H| T + 2(n s H 2 + n T H 2 s )} - 2(D - l)r D H T H s H ST ] 

2r2D _ lH 2+n T /(D-2) H 2+,i s /(D-2) ' 



One finds that the the above curvature invariant quan- 
tity necessarily blows up at (|3 .41 [) . where the inverse of 
the lapse function vanishes. Since the circumference ra- 
dius vanishes at (|3.41[) . these two singularities are both 
sitting at the physical center R — 0. The other curva- 
ture scalar quantities (e.g. TZ a bcdTZ abcd and 1Z) do indeed 
diverge at these central points. 

It is notable to remark that the t = surface is 
completely regular, since the curvature invariants are 
bounded therein. The nonzero U(l) charge smoothenes 



the big-bang singularity which exists in the background 
expanding universe (see Figure[T|). Note also that the null 
surface r = is not singular, as in the case of [ll|. We 
may therefore extend the spacetime across the r D ~ 3 = 
surface, i.e., the negative r region for even dimensions, 
and the r 2 < region for odd dimensions. 

Accordingly, the allowed coordinate regions are 
H^ T Hg S > 0, which is arranged to give 

* > t s (r) with r°- 3 > -Q (3.44) 
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or 

t<t s (r) with r D ~ 3 <-Q. (3.45) 

We are not interested in the latter domain since it turns 
out to be disconnected with the outside region r > by 
the singularity. For this reason, we shall focus on the 
coordinate ranges 

t>t.(r), r°- a >-Q, (3.46) 

Next, it is important to address the causal structure of 
these singularities. We can intuitively expect that both of 
these singularities have the timelike structure, since the 
Misner- Sharp mass is unbounded below as approaching 
these singularities. To establish the timelike signature of 
singularities, it is sufficient to show that there exist an 
infinite number of radial null geodesies stemming from 
and terminating into the singularities (Til [43[ . 

In the vicinity of the singularity r = —Q 1 ^ ^ 3 ^, let us 
assume the following asymptotic form of the radial null 
geodesies, 

r + Q 1/(D - 3) =C (t-ti) 9 , (3.47) 

where q(> 0) and Co are constants. Substituting this 
ansatz into the outgoing radial null condition dr/dt = 

2 



c i +ns /2 = Hs + 2 Qns/2 ^ _ nT (3 48) 

This verifies that infinite geodesies parametrized by the 
departure time t\(< to) emanate from the singularity. 
For the ingoing null geodesies, we can also arrive at the 
conclusion that there exist infinite ingoing null geodesies 
running into the singularity. The proof proceeds similarly 
for the singularity t = t s (r) : it is shown that there exist 
infinite ingoing and outgoing geodesies falling into and 
emerging from the singularity at the time t > to for r < 
and at the time t < for r > 0. It follows that both 
of these two singularities have a timelike structure, i.e., 
these are locally naked singularities. As argued in (23| . 
these timelikc singularities are associated with the U(l) 
charge, not with the background universe. 

In order to verify whether these singularities are visi- 
ble to observer far in the distance from the central region, 
we have to trace the geodesies all the way out to infin- 
ity. We will numerically confirm in later section that the 
singularities t s (r) and r — — Q 1 /( D_3 ) lying in the r < 
region are contained inside the black hole. On the other 
hand, we will see that the singularity t s (r) emerging in 
the r > region is not covered by the event horizon. 

IV. SPACETIME STRUCTURES 

We shall explore the spacetime structure of the "black 
hole candidate" (|2.12|) . In order to conclude that the 



present spacetime is indeed qualified as a black hole, we 
need to appreciate the infinite future of the spacetime to 
assess the locus of the event horizon. Thus, the notion of 
event horizon is of less practical use especially in the time- 
evolving spacetimes. In the dynamical spacetime where 
the null geodesic motions are not solved analytically, it 
is much more advantageous to focus on the trapped sur- 
face, which is locally defined and hence free from such 
a difficulty. When a catastrophic gravitational collapse 
occurs to form a black hole, there arises a region where 
even "outgoing" null rays are dragged back due to strong 
gravity of the black hole. These null rays with negative 
expansion will be absorbed into the hole, characterizing 
the intuitive idea of the black hole as a region of no es- 
cape. For each time slice, this surface defines an apparent 
horizon [3?) as an outermost boundary of the trapped 
region in the asymptotically flat spacetimes. Hayward 
generalized these quasilocal notions to define a class of 
trapping horizons [44| . Trapping horizons are associated 
not only with black holes but also with white holes and 
cosmological ones, suitable for the present context. 

Since the event horizons and trapping horizons are con- 
ceptually different, there is no a priori relationship be- 
tween them (the black hole in stationary spacetime is 
exceptional, for which the event horizon and the trap- 
ping horizon coincide [UHII). Nevertheless, we can rely 
on the Hawking's theorem that all trapped surfaces ap- 
pear inside a black hole under the null energy condition 
provided the spacetime is asymptotically flat with some 
additional technical assumptions (see Proposition 9.2.8 
in [U for the proof). Although the present spacetime is 
not asymptotically flat, taking into the fact that the met- 
ric is well-behaved off of the central singularity at R = 0, 
this criterion is of use for our study. This expectation will 
be backed by numerical analysis of null geodesic motions. 

For simplicity of our argument, we shall focus our at- 
tention only to the 4-dimensional case. We also assume 
to > and Q > 0. We define dimcnsionlcss variables 
t := t/to, f := r/Q and work in the dimcnsionlcss metric 
Q~ 2 ds 2 =: ds 2 = g tlu dx^dx v to simplify our discussion. 
To be specific, we work in the metric 

dl 2 = -T 2 ZdP + S" 1 [df 2 + f 2 (d9 2 + sin 2 6>d0 2 )] , 

(4.1) 

where 

T .= h w = if (4 2) 

' Q ' " [(tf+l)"^(f+l)"s]l/2 ' K ' ' 

with 

n T + n s = 4, (4.3) 

for which the dimcnsionlcss circumference radius is de- 
fined by 

R= [(tf + l)" r (f+l)" s ] 1/4 . (4.4) 

The metric is parametrized by two dimcnsionlcss param- 
eters r and Ut- In what follows, dimcnsionlcss variables 
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associated with the metric (|4.1[) will be denoted with 
tilde. 

The physical meaning of r is naively given as follows. 
The energy density of the scalar field at t = to is p^ cx 
I /tl [Eq. ([313]) ] . While, for the extreme RN black hole 
with the charge Q, the total energy density of the U(l) 



2 (cm) I 



fields evaluated on the horizon is given by k p 
1/Q 2 [Eq. dSHSD]- Thus for the time-dependent black 
hole we can claim that r is related to the ratio of two 
energy densities at the event horizon as 



P 



(cm) 



P 



(4.5) 



A more rigorous treatment shown below [see Eq. (|4.44l) ] 
confirms that the above estimate is indeed true. 

Before going to the detailed analysis, we classify our 
spacetimes into nine types [from Case I-(i) to Case III- 
(iii)]. Classification of Cases I— III is dependent on the 
expansion rate of the universe, i.e., whether the universe 
is decelerating (ny < 2), marginal (ut = 2) or accelerat- 
ing (tit > 2). The parameter r is also important in the 
subclassification for each Class I— III, whereby the global 
structures of the solution change drastically. We sum- 
marize our classification in Table [H which includes some 
important results given later. 



however, that values of expansions themselves are not 
invariant quantities due to the relative normalization of 
and . Only their product 



/+(/_ 



-2R~ 2 (V ^R)(V^R) 



(4.9) 



has a covariant meaning. 

A metric sphere is said to be trapped (untrapped) if 
9+9- > (9+9- < 0), and marginal when 6+9- = 0. 
A trapping horizon is the closure of a hypersurface foli- 
ated by marginal surfaces, i.e., trapping horizons occur 
at 6± = 0. 

It is instructive first to see the case of FLRW universe 
(Q = 0), in which only 9- = has a solution 



fTH 



4to 
TlT 



l-n T /2 



(4.10) 



Using the relations (|2.25[) and (|2.28|) . the trapping hori- 
zon is expressed in terms of the cosmic time t as tth = 
p~ 1 io(t/ia) 1 ~ p . Written in the the circumference radius 
R = ar, this is translated into 



(4.11) 



A. Trapping horizons 

A trapped surface is defined on a closed and orientable 
spacelike surface S. Since the present spacetime admits 
an SO(3)-symmctry, we take 5* as a metric sphere. Let 
us introduce the Ncwman-Penrosc null tetrads as 



n^dx^ = 



'-(-rdt + S _1 df) , 



-(-rdi-s-Mf), 



(4.6) 



'2E 



(d9 + i sin 



with m M being a complex conjugate of m^. These null 
vectors satisfy I^ti^ = — 1 



-m^m M . Z M and are the 



future-directed radial null vectors. Using these tetrads, 
expansions are defined by 



+ := 2m i > 1 m l ' ) V t ,l u 



L := 2?ti (m to 1/) V ai ?v , 



(4.7) 



where the operation of "Re" has been omitted because 
the vanishing rotation is obvious. The present met- 
ric (|4.ip computes to give 



n T f 2 Hs^H^H^ ± T (A?H T Hs - H ST ) 
d ± = ^ ZZZTT^ITiaZITT^Ta ■ ( 4 -°J 



2V2rr 2 H, 



which is the Hubble horizon as desired. This kind of trap- 
ping horizon is associated to the background universe-the 
past trapping horizon in the terminology of Hayward (40|- 
manifesting a situation that even the ingoing null rays are 
pushed outwards due to the cosmic expansion. 

Unfortunately, for Q ^ it is not feasible to obtain the 
orbits of trapping horizons in the explicit analytic form 
t = £th(^) unless ht takes an integer. Reminding that 
the trapping horizons occur where V M i? becomes null, 
we find that it is more helpful to express the positions 
of trapping horizons in terms of the circumference radius 
R. Now, t is expressed by R and f as 



R 4 



(f+ l)™s 



l/n T 



(4.12) 



Inserting this relation into Eq. (|2.12l) . we can write the 
metric in terms of f and R as 



ds 



2 --T 2 



dr f R^ nT (n T + if) 
Jr I ~ n T (f + 1)4/«t 



AR ns / nT dR 



1 2 



n T R(f + l) n s/^T 



%df 2 + R 2 dn 2 



(4.13) 



s 



These expansions characterize the rate of divergence of 
outgoing and ingoing light rays. It should be remarked, 



Setting R = constant, one obtains the induced metric of 
the constant R surfaces, 
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TABLE I: Classification of the spacetime structures of the solution (|4,l|l . The parameter tit corresponds to the steepness 
parameter of the dilaton potential, which is related to the expansion exponent of the FLRW universe a oc i p as Eq. (|2.28|) . while 
r := to/Q measures the "nonextremality" of the black hole. This classification is based on the spacetime structures discussed 
in the text. In the last column, "BH" and "NS" represent whether the spacetime describes a regular black hole or a naked 
singularity without an event horizon. 



Case 


TIT 




number of 


causal structure 


T 


horizons 


infinity 


near-horizon 


global structure 


extreme RN 


TIT = 


any 


1 


Minkowski 


Fig. [2] (left) 


Fig- M (right) : BH 


I 

(decelerating universe: p < 1) 




< n T < 4/3 






Fig. 


Hi) 








n T = 4/3 


any 


2 


Fig. 


[1121 


Fig. ETA) 


Fig. Eta) : BH 


(iii) 


4/3 < tit < 2 






Fig. 








II 

(Milne universe: p = 1) 


0) 




T > 1 


2 






Fig. ETA) 


Fig. Eta) : BH 


(ii) 


TIT = 2 


T = 1 


1 


Fig. EM 


Fig. [HE) 


Fig. [3» : NS 


(iii) 




r < 1 






Fig. lie) 


III 

(accelerating universe: p > 1) 


(0 




t > r cr 


3 






Fig. MP) 


Fig. Eld) : BH 


(ii) 


2 < n T < 4 


T = T cr 


2 


Fig. [115) 


Fig.[5TE) 


Fig. I7[e) : BH 


(iii) 




T < Tcr 


1 






Fig. IMC) 


Fig. Etc) :NS 


RNdS 


TIT = 4 


any 


3 


dS 


n/a 


Fig. MP) : BH 



d§l 



R 2 r 2 



J R 4 /" T (n T + 4f) 
n T (l + f) 4 /«T 



i? 4 



df 2 , 
(4.14) 



where we have discarded the angular parts which are ir- 
relevant to the discussion of trapping horizon. 

Figure [3] depicts the R — f diagram, from which we 
can read off the signature of R = constant surfaces. 
From Eq. (|4.9[) . the region where the R = constant 
surface is timelike (spacelike) consists of the untrapped 
(trapped) surfaces. The thick lines in Figure [3] desig- 
nate the loci of trapping horizons at which the terms 
in square-bracket of Eq. (|4.14[) vanish, dsj^ = 0. The 
spacetime has at most three types of trapping horizons, 
Rz(r) < i?i,_(f) < Ri t +(r), where R,2(f) satisfies 



^ 2 (r) 4/ " T (nT+4f) _ 
n T (l+f) 4 / n T 



+ R 2 {ff = {), (4.15) 



whereas i?i._(r) and i?i j+ (f) are two positive roots (if 
they exist) of the equation 



£i(f) 4 /* 



(n T + 4f) 



n T {l + r) 4 / n T 



R 1 (f) 2 = 0. (4.16) 



Observe that the definition of Ri(f) and R2(r) differs 
from that in the previous paper [ll|, where Ri(r) and 
R2(f) were defined by the surfaces of (9_ =0 and 9 + = 0, 
respectively. In the f < region, the surfaces R± (f) and 
R-2(f) defined by Eqs. (|4.16[) and (|4.15p arc not always 
coincident with the 6L = and 9 + = surfaces. Still, it 
is found that in the f > region = holds at Ri(f) 
and 9 + = holds at R.2{r). 



In Figure [3] we display seven plots of trapping horizons 
according to the classification shown in Table |U [We 
have dropped Cases I- (ii) and I- (iii) because these plots 
are the same as that for Case I-(i).] Our previous pa- 
per [ll| discussed the spacetime with tit = 1 belong- 
ing to Case I-(i), where two trapping horizons i?i > _(r) 
and R2{r) arise. Case II-(i) behaves similarly as Case I. 
The universe in Case II undergoes a marginal accelera- 
tion a oc i. Our main interest is toward Case III, where 



the background FLRW universe is accelerating. Figure [3] 
illustrates that the region R > i?i, + (f) with r > is 
composed of trapped surfaces, in accordance with the 
acceleration of the background FLRW universe. 

Let us uncover the properties of these trapping hori- 
zons. Looking at Eq. (|4.15j) . one finds that the trapping 
horizon R2 (f) approaches to the constant R2 (f) — > y/r in 
the limit of f — > 00. Noting the inequality dR,2(r)/dr > 
for r > 0, i?2 (f) < ^/t is concluded in the f > region. In 
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Case I-(i) 




Case II-(i) 
R + 



.0 -05 0.0 0.5 1.0 U 2.0 



7?_ timelike 



Case Il-(ii) 



sp.icelikc 



Case Il-(iii) 



Space like 



- 1.0 - 0.5 0.0 0.5 1.0 1.5 2.0 



■1.0 -0.5 0.0 0.5 1.0 1.5 2.0 



Case III-(i) 







R c 






timelike \ 












spacclike 




timelike 



Case Ill-(ii) 



R c - R + 



spacclike 



Case Ill-(iii) 



1.0 -0.5 0.0 0.5 1.0 



1.0 - 5 0.0 5 1.0 1.5 2.0 



spacelike 



R- timelike 



1.0 - 0.5 0.0 0.5 1.0 1.5 2.0 



FIG. 3: Typical plots of circumference radius R (vertical axis) against f (horizontal axis) for Case I-(i) {tit = 1 with r = 1), 
Case II-(i) (n T = 2 with r = 3/2), Case Il-(ii) (n T = 2 with r = 1), Case Il-(iii) (n T = 2 with r = 1/2), Case III-(i) 
(tit = 3 with r = 3), Case Ill-(ii) (tit = 3 with r = r cr ) and Case Ill-(iii) (tit = 3 with r = 1). The curve of Case I-(i) 
is the representative for Case I. The thick lines denote trapping horizons, across which the R = constant surface changes the 
signature. The region "spacelike" (and respectively "timelike") stands for the domain where the R = constant surfaces are 
spacelike (timelike). In Case III-(i), the bounded region tends to shrink as r decreases and eventually degenerate into a point 
shown in Case Ill-(ii). 



Case I {n T < 2), dR/df < (>)0 holds for R 2 < (>)R 2 ,< 
in the region f < 0, where 



R 



2,< 



2t 



2-n T 



(4.17) 



While, for ut > 2, cLffe/df < holds in the r < region. 
Thus, in every case the trapping horizon R2(f) develops 
and diverges as f — >• — ny/4. 

It is easy to find that i?i i +(f ) is confluent with (f) 
at the constant radius R = Ri t >, where 



Ri.> ■= 



2t 



(4.18) 



at which Ri,±(r) has a vertical tangent. Hence, the 
largest trapping horizon i?i i+ (f) exists if and only if 
tit > 2. For tit < 2 and tit = 2 with t > 1, 
the trapping horizon i?x,-(f) diverges as r — > —tlt/^. 
For n T > 2, dfli _/df > (<)0 and d^i,+/df < (>)0 
hold in the region f > (<)0. Denoting i?_ := i?2(0), 
R + := Ri-(0) and R c := Ri,+ (0), it then follows that 
R+ < Ri-(f) < Ri,+ (f) < R c for Case III-(i). 

The signature of trapping horizons is an important is- 
sue to be discussed. Differentiating Eq. (|4.15p . we ob- 
tain the expression of AR 2 / df along the trapping horizon. 



Substituting this into the metric (|4.13|) . we can find the 
induced metric on the trapping horizon. After a simple 
but tedious algebra, we obtain 



4(1 



nsr 2 R 2 



n T R\ 



2R 2 2 (t - Rl) 



dR 2 2 . (4.19) 

{tit + vr) 1 

It immediately follows that the cases of tit > 2 and tit < 
2 with R < R 2y< assure that the trapping horizon R 2 (f) 
is always timelike d.s^ H < 0. We can numerically check 
that i?2,- is timelike also for R > i?2,<- 

Similarly, along the trapping horizon R±, we have 



4(1 



nsf R\ 



n T R\-2Rl(T + Rl) 



4n s f 2 (T + Rl) 2 



(tit + 4f) 



P\2 



dR\ . (4.20) 



Inspecting Eq. (|4.18[) . the trapping horizon R\ t + is nec- 
essarily timclike. On the other hand, to see the signature 
of R\-, numerical calculation is necessary. Nevertheless, 
we can see the universal behavior around f = 0, where 
the last term in the square bracket of Eq. (|4.20p is neg- 
ligible. Taking into account Eq. (|4.18p . it is found that 
Ri - is spacelike near f = 0. 
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An analytic estimation is possible for the background 
FLRW universe Q = 0, for which Eq. (|4.20p reduces to 



rU 2 

as TH 



4 - 3n T 
ns 



dRf, 



(4.21) 



implying that the signature of trapping horizon is time- 
like for p > 1/2 (tit > 4/3), spacelike for p < 1/2 
(n T < 4/3) and null for p = 1/2 (n T = 4/3). This ac- 
counts for the subclassification of Case I shown in Fig. [1] 
The behaviors of trapping horizon Ri — in the asymptotic 
region are analogous to those in the FLRW case. 



where T> ^ is a derivative operator of the near-horizon 
metric (|4.23j) . Hence, by transforming to the new coor- 
dinates 



where 



T = ±ln|i| 



R 4^3+4/r lrd ^ 



R= (l + tf) nT ^, (4.26) 



f(R) :=t 2 (k i/nT -l) - R 



(4.27) 



B. Event horizons 

In order to identify the locus of event horizon, we follow 
the strategy laid out in our previous paper [ll| . We shall 
first discuss the geometry of the near-horizon metric and 
demonstrate that the null surfaces of "event-horizon can- 
didate" in our dynamical spacetime (|4.ip are described 
by the Killing horizons. Our next tactics is to confirm 
these local horizons arc indeed qualified as true horizons 
in the original spacetime by computing the null geodesic 
motions numerically. These logical steps will guide us in 
the right direction for obtaining the global causal struc- 
ture. 



1. Near-horizon geometry 

The analysis in Sec. [TT| reveals that the limit r — > 
and t being finite corresponds to the "throat" geometry. 
This fact leads us to speculate that the event horizons 
in the present spacetime, if exist, correspond to the null 
surfaces at r — > with i —> ±oo, i.e., the infinite rcdshift 
and blueshift surfaces "joined" at the throat. To discuss 
these null surfaces in more detail, it is useful to take the 
near-horizon limit, defined by 



t 



t 



cr . 



0. 



(4.22) 



After the rescalings, the metric is free from e, and we can 
obtain the near-horizon geometry 



d.s 



NH 



-T^ 2 (l+tf)-^/ 2 dP 

+ f' 2 (1 + tf) riT/2 (df 2 + f 2 dn 2 2 ) . (4.23) 



As a direct corollary of the scaling limit (|4.22[) , the near- 
horizon metric (|4.23j) admits a Killing vector 



(4.24) 



By a straightforward calculation, we find that the above 
Killing field is hypersurf ace- orthogonal 



i{yP v £ p ] = . 



(4.25) 



we can bring the near-horizon metric (|4.23[) into a man- 
ifestly static form, 



d.s 



7,2 



f(R) 

R 2 



dT 2 



l6T 2 R S/n T 

n 2 T f(R) 



dR 2 



R 2 dn 2 . (4.28) 



In this coordinate system, the Killing vector (|4.24[) is 
expressed as £ M = (d/dTY . It is obvious that the above 
metric (|4.28p describes a static spacetime with Killing 
horizons at f(R) = where the Killing vector £ M becomes 
null. It deserves to note that terms in square bracket in 
Eq. (|4.14[) simplify to f(R) in the near horizon limit. 
This means that the trapping horizons tend to be null as 
approaching f — > 0, in accordance with the results in the 
previous subsection. 

Since t and f are not gauge-invariant quantities, one 
may be worried about whether the scaling limit (|4.22j) is 
indeed well-defined. This prescription of "zooming-up" 
can be justified as follows. In the scaling limit (|4.22p . the 
dilaton and the U(l) fields reduce to 



K<f> 



2n s 
riT 



InR, 



and 



K F^ = 

TR 



W2TTT 



4V27TT 



(4.29) 



(4.30) 



We can confirm that the spacetime (|4.28|) with (|4.29p 
and (|4.30p still satisfies the field equations of Einstein- 
"Maxwell" -dilaton gravity (|2.3p ~ (|2.5p . This warrants 
that the near-horizon limit (|4.22p is not the singular pro- 
cedure since it maps the original metric (|2. 12|) onto an- 
other one (|4.23p in the same theory. Note that the case 
n-r = 4 does not yield the near-horizon metric, since the 
metric (|4.23p is identical to the original one (|4.ip and 
simply giving rise to a static chart of the RNdS black 
hole with M = Q. 

Let us devote some space to discuss the near-horizon 
geometry in more detail. Although the near-horizon met- 
ric (|4.23p is neither asymptotically flat nor asymptoti- 
cally AdS, the global causal structures reduce to those of 
known solutions. Such unusual asymptotic structures of 
black holes are commonly found when the dilaton field 
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is present (46l - |48j . The near- horizon metric (|4.23j) ad- 
mits horizons at f(R) = 0, which is rewritten into the 
equation 

F ± (R) = r, with F ± (R) := — g^- . (4.31) 

This equation gives the candidates for horizon radius, 
J?_, and R c , which were found in the limit of r — > 
of the trapping horizons in the previous subsection. 
The number of roots for Eq. (|4.31|) [or f(R) = 0] can 
be seen visually as the intersection points of F±(R) and 
t = constant (see Figure [4} . 

Irrespective of the values of r and tit, behavior of F— 
is universal: it starts from zero and diverges as R — > 1. 
Then, F_ (R) = r has always a solution R = ii_ less than 
unity. Whereas, the behavior of F + depends on the value 
of tit, i.e., rir < 2, tit = 2 or tit > 2. As is clear from 
Figure QJ for tit < 2 [Case I in Table |T] and tit = 2 with 



t > 1 [Case II-(i)], we find two horizons (i?_ and i?+). 
When tit = 2 and r < 1 [Cases H-(ii) and (iii)], we have 
only one root {R-)- In the case of tit > 2, where the 
background universe is accelerating, F + has a minimum 
at 

(\ n T /4 
^h) • (432) 

Hence, for 

r > r cr := \n n T T/ \n T - 2) 1 -"-/ 2 . (4.33) 

the equation (|4.31[) has three distinct roots, R_ < R + < 
R c [Case III-(i)]. When t = r cr , two horizons i?+ and 
R c become degenerate [Case Ill-(ii)]. If t < r cr , we have 
only one root [Case II- (iii)]. These results are shown in 
Table D 




Case I Case II Case III 



FIG. 4: Plots of F+ and F- for Case I (left), Case II (middle) and Case III (right) (see Table U for classification). The equality 
f(R) = then possesses two distinct solutions for Case I and Case II-(i), a single solution for Cases Il-(ii), II- (iii) and III- (iii) , 
three distinct solutions for Case III-(i) and two solutions with larger one being degenerate for Case Ill-(ii). 



The spacetime structure of the near-horizon metric for 
each case I— III is obtained by simple evaluation of radial 
null geodesies (see e.g., HH). We will not devote much 
time to discussing this issue since our primary interest 
in this subsection is exclusively to the horizon structure. 
We obtain 5 types of conformal diagrams, as shown in 
Figure [5] (A)-(E). The figure (A) corresponds to Cases I 
and II-(i), (B) to Case H-(ii), (C) to Case II- (iii) and III- 
(iii), (D) to Case III-(i) and (E) to Case Ill-(ii). These 
classifications are also included in Table |U Despite the 
nontrivial asymptotic structures (|4.28p . these conformal 
diagrams are identical to the familiar solutions. The fig- 
ure (A) is the same as the RN-AdS, (B) as the vanishing 
mass RN spacetime with a hyperbolic horizon, (C) as the 
negative mass SdS, (D) as the nonextremal RNdS and (E) 
as the degenerate RNdS, respectively. The surfaces 



and i?_ are the black hole event horizon and the white 
hole horizon for the near-horizon geometry and also for 
the original metric. It should be emphasized, however, 
that the largest horizon R c does not approximate our 
original metric, since the surface R c corresponds to the 
limit t — > and f — > oo with tf being finite in the near- 
horizon geometry (the green line in Figure [S]) . It inci- 
dentally arises due to the symmetry of the near-horizon 
metric not shared by the original metric. Consequently, 
the null surface R c is irrelevant to the cosmological hori- 
zon in the original metric (|4.1[) , although it is a cosmo- 
logical horizon for the near-horizon metric. It is worth 
commenting also that the case tit — ns — 2 with r = 1 
is special, for which there exists "internal null infinity" 
lf n specified by t — > oo and f — > with R= (tr + l) 1 / 2 
being infinite. This surface has an outgoing null struc- 
ture. Only the ingoing null geodesies can approach lf nl 



16 



hence it should be distinguished by X^~ ut (t oo and 
f — > oo). The Carter- Penrose diagram in this case is 
given in Figure [5JB). 

To conclude, the neighborhoods of the null surfaces R+ 
and R- in the original metric (|4.1[) are locally isomet- 
ric to the corresponding null surfaces-these are indeed 
the Killing horizons-in the near-horizon metric (|4.23|) or 
(|4.28j) . In Figurc[5l only the shaded regions in the vicinity 
of the horizon approximate our original metric (|4.ip . The 
surface if exists, deserves a black- hole event horizon 
for the near-horizon metric and it turns out later, this 
surface corresponds to a black hole horizon for the origi- 
nal metric as well. 

Note also that the vector field ^ (|4.24l) solves the 
Killing equation in the original spacetime only at the 
horizon. Since the outside region is highly dynamical, it 
comes as a surprise for us that the unexpected Killing 
symmetry arises. This is achieved due to the fact that 
the electromagnetic repulsive force and the cosmic expan- 
sion are compensated by the gravitational and dilatonic 
attraction. These Killing horizons R± are in general non- 
degenerate in the sense that the surface gravities k± are 



nonvanishing. From the formula 

24=T(V^)(V /i r), (4-34) 
one obtains the surface gravities associated with £ M as 



k± = It 



n T R 2 ± i/nT 
2t 



(4.35) 



These are nonvanishing unless r = r cr . The surface grav- 
ities are constant over the horizon, as expected. 



2. Null geodesies 

The argument given in previous subsection strongly 
implies that i?+ and i?_ are the future and past event 
horizons. In order to conclude this in a more rigorous 
manner, we need to discuss geodesic motions. Since the 
present spacetime admits a spherical symmetry, it is suf- 
ficient to focus on the radial null geodesies to determine 
the causal structure. 



Using the relation for the radial null vectors, 



rt = ±JH!£ r H s 
the radial null geodesic equations simplify to 
tH S t 



t± 



0, 



■■ n T {H T H s J ~2 _ n 
2tH t 



(4.36) 

(4.37) 
(4.38) 



Here and hereafter, the plus (minus) sign refers to as the 
outgoing (ingoing) null geodesies, and the dot means the 
differentiation with respect to an affinc parameter A. 

Following [ll|, [Ho| the asymptotic geodesic behavior 
around the surface R + (corresponding to t — > oo and 
r — >• 0) can be evaluated as 



i +) (A-A + ) 1 ^+, f=4 +) (A-A + ) 



(4.39) 



where is constant, A+ corresponds to the arrival time 



at R + and k + was defined by Eq. (|4.35p . This equation 
not only means that the horizon can be reached within 
a finite affine time but also that r and l/t are not the 



smooth functions of A at A + (unless k + 1 happens to be 
integral) . 

Let us move on to the numerical analysis. Exploiting 
the affine transformation, if necessary, we can impose 
without loss of any generality A — at the starting point 
of the geodesies and i(0) = ±1. Then the geodesies are 
distinguished by their initial point [t(0),r(0)]. 

Our primary concern here is Case III, where the back- 
ground universe is accelerating. We will give a detailed 
description exclusively of ut = 3 with r = 3 [Case 
III-(i)]. We first develop geodesies for the outside re- 
gion r > by setting the initial radial coordinate to 
be f(0) = 1. It then follows that the r = 1 sur- 
face is divided by three trapping horizons into four re- 
gions (see Figure Take the representative spacetime 
points pi (I = 1,...4) specified by (f/(0),r(0) = 1) with 
ti(0) > t 2 (0) > t 3 (0) > £4(0) in such a way that pi is 
contained in each partitioned (un)trapped region. We 
shall refer to the geodesies starting from the initial point 
pi as "Class-/." For f < 0, we take f(0) = —1/10 and re- 
peat the identical procedure with redefining "Class-/" as 
i/i(0) < t 3 (0) < t 2 (0) < £i(0). Figure [5] presents our nu- 
merical result in the accelerating universe (tit = T = 3) 
and the deduced Carter-Penrose diagrams are shown in 
Figure Eld). 



J r 



Let us begin by the analysis for f > 0. Among 
the Class- 1 geodesies of the future-directed ingoing null, 



geodesies with a sufficiently large t\ tend to diverge as 
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(A) Cases I and II-(i) (B) Case Il-(ii) (C) Cases Il-(iii) and Ill-(iii) 




(D) Case III-(i) (E) Case Ill-(ii) 




FIG. 5: Conformal diagrams of the maximally extended near-horizon geometries. The global structures are the same as those 
for (A) the nonextremal RN-AdS, (B) the zero-mass RN spacetime with a "hyperbolic horizon," (C) the negative mass SdS 
(or the overcharged RNdS), (D) the nonextremal RNdS and (E) the extremal RNdS. The thick dotted-dashed lines denote the 
central singularity at t = —1/r. The shaded regions approximate our original solution (14.111 . 



the affine time evolves, while geodesies with ii(0) being 
close to the trapping horizon converge to i?+ with an in- 
finite rcdshift. This implies the existence of a cosmolog- 
ical horizon [the 45° line denoted by CH in Figure [7(d)], 
which is a characteristic feature in the accelerating uni- 
verse as seen for the p > 1 case in Figure [1] Other 
classes of geodesies unavoidably approach i? + within a 
finite affine parameter. 

Behaviors of future-directed outgoing null geodesies of 
Class-3 imply the existence a critical null curve [the black 



dotted 45° line in Figure [7(d)] above which the geodesies 
can reach infinity and below which the geodesies in- 
evitably plunge into the singularity t = t s (f). 

Motions of the past-directed ingoing null geodesies are 
universal: they necessarily get to the null surface i?_ 
with undergoing an infinite blueshift. 

The past-directed outgoing null geodesies also exhibit 
a universal behavior: they inevitably reach the timelike 
singularity t s (r) [the gray dashed-dotted line R = in 
Figure 17(d)] within a finite affine parameter. 



Assembling these results, we are able to draw a con- 
formal diagram of Figure [7(d). 3 The future infinity con- 



3 We traced several geodesies starting from different initial points 
and checked consistency. 



sists of a spacelike surface, consistent with the behavior 
of future-directed ingoing null geodesies of Class- 1, and 
with the FLRW limit r — > oo as described in Figure QJV). 
The timelike singularity t s (f) = —1/r exists up to t = 0, 
which is joined with X + . The cosmological horizon de- 
velops, whose orbit t = tcn(^) can be traced numer- 
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FIG. 6: Behaviors of radial null geodesies for the accelerating universe (nr = 3 with r = 3) in the outside region (left) and the 
inside region (right). The figures show the circumference radius versus the affine parameter. The red, blue and purple lines 
stand for R- ~ 0.82, 7?+ ~ 1.56 and R c ~ 4.03, respectively. The diagrams in the top, upper middle, lower middle, bottom 
rows correspond to future-directed ingoing null geodesies, future-directed outgoing null geodesies, past-directed ingoing null 
geodesies and past-directed outgoing null geodesies. We have shown the null geodesies for the initial values f(0) = 1 (left) and 
f(0) = -1/10 (right). 



ically. The green-dotted lines denote the R = constant 
surfaces, which change signature across the trapping hori- 
zons. The R = constant surfaces can be written in com- 
parison with Figure [3J and the trapping horizons can be 
drawn from Eq. (|4.19j) and (|4.20|> . The Carter-Penrose 
diagram explicitly shows that the null surface R+ de- 
serves to be a black hole horizon. The null surface i?_ , on 
the other hand, may not suitable to be called a white hole 
horizon in a rigorous sense since the spacctime admits no 
past infinity in the outside region. Nevertheless, we shall 
stick to call it a "white hole horizon," from the spirit of 
a one way membrane as a "region of no-entrance." 

The global structure of the solution in the inside region 
(r < 0) can be deduced analogously. It is notable that 
the inside region admits a past null infinity I~ . At first 
glance, it seems peculiar because the background uni- 
verse is expanding (at X~ the circumference radius R is 
already infinite there and the universe is not able to ex- 
pand further!). We recognize that the spacetime indeed 
approaches to the expanding FLRW universe only in the 
limit of r ->• +oo with t/t > [see Eq. (|2.25j) ]. while 
this is not true in the limit t/to — oo with — 1 < r < 0, 
in which case the background geometry looks like a con- 



tracting universe. 

Following the above treatment, one can draw the 
Carter-Penrose diagrams for other cases. We do not dis- 
play those behaviors of null geodesies corresponding to 
Figure[5]in order to retain the compactness of the present 
paper, but they are easy to obtain. (Without geodesic 
analysis one can intuitively deduce the causal structures 
simply by inspecting the R — f diagram in Figure [3j . 
We have arrived at 5-types of global conformal diagrams 
shown in Figure [7] (a)-(e). Let us summarize the results. 

Figure [7|a) corresponding to Cases I and II-(i) was 
given in our previous paper [ll| . The spacetime has a reg- 
ular, nondegenerate event horizon and approaches to 
the decelerating FLRW universe. In the picture shown 
in Fig [7|a), we sketch the trapping horizons for the 
riT = t = 1 case, where the spacelike trapping horizon 
i?i develops in the outside domain. Its signature depends 
sensitively on n-r and r. For example, it becomes asymp- 
totically null in the radiation-dominated case (ut = 4/3). 
Though the signature of trapping horizons is modified de- 
pending on the parameters, the causal structure remains 
unaltered. 
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(a) Cases I and II-(i) (b) Cases II (ii)-(iii) (c) Case Ill-(iii) 




(d) Case III-(i) 




(e) Case Ill-(ii) 




FIG. 7: Conformal diagrams of the present spacetime (|2.12[) for Cases I— III. The green dotted curves stand for the R — constant 
surfaces. The red thick curves represent the trapping horizons (TH) and blue lines denote the cosmological horizons (CH). The 
gray and black dotted-dashed curves correspond to the singularities of r — —Q and t — t s (r), respectively. Figure (a), where 
the universe is decelerating, has been argued in the previous paper [llj]. For Figure (b), the universe undergoes a marginal 
acceleration a oc i . It admits an internal null infinity 1^ at which r = and t — > oo with R — s- co. Figures (c)-(e) correspond 
to the accelerating universe. Figure (d) shows that there appear three horizons: the black hole horizon R+, the "white hole 
horizon" i?_ and the cosmological horizon denoted by CH (not R c ). In Figure (e), the horizon is degenerate, while in Figure 
(c) the spacetime is overcharged to give a globally naked singularity. Figures (a), (d) and (e) are extensible across the surface 
R~ with a contracting universe in a continuous but nonanalytic manner. 



In Figure [Tfe) [Case Ill-(ii)], the spacetime has a de- 
generate event horizon at R + = R c and asymptotically 
tends to the accelerating FLRW universe. This is a spe- 
cial situation of Case III- (i) . 

On the other hand, the spacetime fails to have a black 
hole horizon in FiguresJTJb) [Cases H-(ii) and Il-(iii)] and 
Etc) [Case Ill-(iii)] (see Tabled]). These cases correspond 
to spacetimes with naked singularities. In Figure [T^b) 
where tit = n$ = 2, one may recognize that there exists 
the "internal" null infinity lf n . where f — > and t — > oo 
with R ~ {If) 1 / 2 — > co. We can show the existence of 
Xf n as follows. Setting tit = ns and assuming f ~ 



with t — > oo, one can find the asymptotic solutions of 
the radial ingoing null geodesies (|4.37p and (|4.38[) around 
A — > oo as 

r = \- T/{1 - T \ with t = CiA 1/(1 ~ T) , (4.40) 
for t < 1, and 

r = E^), with t = X, (4.41) 
A 

for t = 1, where G\ and C2 are positive constants and 
W is Lambert's TV-function satisfying W{x)e w ^ = x. 
These asymptotic geodesies both satisfy tf — J- 00 as A — > 
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FIG. 8: Relation between r and the ratio of two energy den- 
sities at the horizon in the double logarithmic plot. This 
illustrates the almost linear relationship between r 2 and 

IP '\H- 



oo. Whereas, the outgoing radial null geodesies do not 
decrease f. This proves the existence of I m , at which 
only the ingoing null geodesies can arrive. 

We have revealed that our solution (|4.1[) specified 
by three continuous parameters (Q,tit,t) describes a 
charged black hole in the FLRW universe in some pa- 
rameter range. It is therefore worthwhile to discuss the 
physical meaning of these parameters. It is clear that 
the constant Q is a U(l) charge, whilst ttt is the steep- 
ness of the dilaton potential (if tit happens to be natural 
number, it counts the number of time-dependent branes). 

If riT = 4, the parameter t$ has a clear meaning as 
a curvature radius of the dc Sitter spacetime. However, 
the meaning of r (or to) is less obvious for tit ^ 4. In 
order to clarify the meaning, we shall consider the energy 
densities of two fields. The values of p^ and p( om ) on 
the event horizon i? + are given by 

»t 



1 1 4Q 2 i? 2 



-8/n T 



which leads to 

iM 2r 2 



P 



H nxiriT + 2) 



i?7 4 n T + n s R 



8/n T 



(4.43) 



(4.44) 



The horizon radius R + is given in terms of r via 
Eq. (|4.3ip . As a result, we find that the ratio of two 
energy densities at the horizon is related to r, although 
the relation is complicated. In the limit of r — > oo, how- 
ever, we find a simple relationship 



(cm) 



P 



(*) 



8 



tit(tit + 2) 



(4.45) 



Figure [5] depicts the dependence of p^ em '/p^' against t. 
It follows that t 2 has a one to one correspondence to the 
ratio of two energy densities evaluated at the horizon i? + . 



We can also give an afterthought implication to r as 
follows. From Eq. (|4.31[) , one obtains a relationship 



R 2 



(4.46) 



This equation implies that the parameter r measures the 
"noncxtremality" of the Killing horizon. An extremal 
limit is taken r — ¥ oo where i?+ = i?_ and r — > r cr where 
R + = R c . Inferring from the supersymmetric case, the 
nonextremality indicates that the mass is strictly larger 
than the (central) charge. Hence, we deduce that some 
combination of the above parameters will give a well- 
defined mass of a black hole larger than the charge Q. 
Considering that the black hole horizon does not occur 
for small r, this point of view is supported. If a well- 
behaved mass definition is found, we arc able to discuss 
the thermodynamics of a black hole. This expected pos- 
sibility is now under investigation. 



V. CONCLUDING REMARKS 

We have presented a family of solutions in a D- 
dimensional Einstein- "Maxwell" -dilaton system with a 
Liouvillc potential. In the single mass case, the solution 
intcrlies the extremal black hole and the FLRW universe 
with a power-law expansion. This solution reduces in the 
special cases to the extremal RN black hole (tit — 0), the 
nonextrcmal RNdS black hole (ns = 0) and the solution 
of derived from the intersecting M-brancs (D = 4 
and tit = 1). The present spacetime is characterized by 
three parameters: the steepness of the dilaton potential 
tit, the U(l) charge Q and the "nonextremality" (or the 
ratio of two energy densities at the horizon) r. The pa- 
rameter tit controls the expanding power p of the back- 
ground FLRW universe. For any choice of parameters, 
the system is shown to obey the weak energy condition. 

The primary aim of this paper is to obtain global struc- 
tures of the solution. The spacetime can be grouped 
into the nine cases [Cases I-(i) — III- (iii)] summarized in 
Table HI according to the parameter values. Our conse- 
quence is visually captured in Figure[71 where the 5-types 
of global structure are obtained. In the case where the 
background universe is accelerating (tit > 2) with large 
nonextremality parameter (r > r cr ), the spacetime in- 
deed describes a charged black hole in the FLRW universe 
undergoing an accelerated expansion [Case III-(i)]. We 
have also clarified the global structures for the marginally 
accelerating and decelerating cases. This has been an 
open issue in 
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The present solution displays various interesting fea- 
tures. At first glance, the solution appears to have a 
degenerate horizon. We have shown that this is not the 
case. Taking the near-horizon limit (|4.22p . the horizon is 
not degenerate except for r = oo and r = r cr . Amazingly, 
the event horizon of a black hole constitutes a Killing 
horizon. This means that the matter fields fail to ac- 
crete into the black hole, whereby the area of the black 
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hole keeps constant. Hence, the present spacetime is not 
suitable for describing a growing black hole in FLRW 
universe. Instead, the solution preserves the analogue of 
the equilibrium state despite the time-dependence of the 
metric. The situation is closely analogous to the super- 
symmetric state. 

The black hole thermodynamics in the expanding uni- 
verse is an interesting future work to be discussed. Since 
the event horizon is generated by a (partial) Killing field, 
the black hole has a nonvanishing surface gravity. How- 
ever, it is far from obvious to which observers the physical 
temperature is assigned, since the surface gravity is sen- 
sitive to the normalization of a Killing field. This may 
be rephrased as what is the mass of the black hole. We 
hope to visit this issue in a separated paper. 

Another interesting issue to be explored in the time- 
dependent black hole spacetime is a black hole merger. 
Adding multiple point sources and reversing time back- 
wards, the present metric turns out to give the ana- 
lytic description of black hole coalescence, similar to 



the Kastor-Traschen solution [2^, [2(|. In the single- 
mass Kastor-Traschen spacetime (RNdS), an analogue of 
"thermal equilibrium" is realized due to the fact that the 
temperature of the event horizon and that of the cosmo- 
logical horizon become the same. If this viewpoint con- 
tinues to be valid in the present spacetime, what plays the 
role of a box containing the black hole which is thermal 
equilibrium with a bath of radiation? These are interest- 
ing questions to be explored. 
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